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First Observation of Diboson Production in Hadronic Final State at the 
Tevatron 

J. Pursley, on behalf of the CDF Collaboration 

Department of Physics, University of Wisconsin-Madison, Madison, Wl, 53706, USA 

We present the first observation in hadronic collisions of the electroweak production of vector boson pairs 
(VV; V = W,Z) where one boson decays to a hadronic final state. The data correspond to 3.5 fb" 1 of 
integrated luminosity collected by the CDF II detector in pp collisions at a center-of-mass energy of 1.96 TeV. 
Event selection requires two identified jets and large transverse momentum imbalance. The analysis employs 
several novel techniques to suppress multijet background and reduce systematic uncertainties. We observe 
1516 ± 239(stat) ±144(syst) diboson candidate events and measure a cross section of cr(pp — » VV + X) = 
18.0 ± 2.8(stat) ± 2.4(syst) ± l.l(lumi) pb, in agreement with standard model (SM) expectations. 



1. Introduction 

The production of heavy gauge boson pairs (WW , 
WZ, and ZZ) has been observed in fully leptonic final 
states at the Fcrmilab Tevatron collider [TJ [2J. Evi- 
dence for diboson decays into an tvqq final state, with 
I = e,[i,T and q = u,d, s, c, b, was recently presented 
by the DO collaboration [3]. The analysis presented 
here is the first conclusive observation at a hadron 
collider of diboson production with one boson decay- 
ing into leptons and the other into hadrons jj. 

Measurements of the diboson production cross sec- 
tions provide tests of the self-interactions of the gauge 
bosons, and deviations from SM predictions could in- 
dicate new physics [5]. Diboson production involving 
hadronic decays is also topologically similar to asso- 
ciated Higgs boson production, pp — * VH + X, when 
the Higgs boson decays to bb, the most promising dis- 
covery channel for a low-mass Higgs boson. Thus the 
analysis techniques described here will also be useful 
for Higgs boson searches. 



2. Experimental Apparatus 

The CDF II detector is described in detail in 
Ref. [6]. Surrounding the beam pipe is a tracking 
system consisting of a silicon microstrip detector, a 
cylindrical drift chamber, and a solenoid producing a 
1.4T magnetic field along the beam axis. The central 
and forward calorimeters surround the tracking sys- 
tem with a projective tower geometry. The calorime- 
ters are composed of inner electromagnetic and outer 
hadronic sections that consist of lead-scintillator and 
iron-scintillator respectively. A system of Cerenkov 
counters, located around the beam pipe and inside 
the forward calorimeters, measures the number of in- 
elastic pp collisions per bunch crossing and thus the 
luminosity [7J. 

The geometry of the detector is characterized by 
the azimuthal angle <fi and the pseudorapidity r/ = 
— ln[tan(0/2)], where 9 is the polar angle measured 



from the proton beam direction. The transverse en- 
ergy Et — E sin 6, where E is the energy in the 
calorimeter towers associated with a cluster of energy 
deposition. Transverse momentum p^ is the track 
momentum component transverse to the beam-line. 
The missing transverse energy vector $ T is defined as 
— E l T h T , where h T is the transverse component 
of the unit vector pointing from the interaction point 
to the energy deposition in calorimeter tower i. This 
is corrected for the pr of muons, which do not deposit 
all of their energy in the calorimeter, and tracks which 
point to uninstrumented regions in the calorimeter. 
The missing transverse energy Ifi T is defined as \]fi T \. 
Strongly interacting partons undergo fragmentation 
that results in highly collimated jets of hadronic parti- 
cles. Jet candidates are reconstructed in the calorime- 
ter using the JETCLU cone algorithm [5] with a cone 
radius of 0.4 in (r),(f)) space. The energy measured by 
the calorimeter must be corrected to improve the es- 
timated energy [9]. The effects corrected for include 
the non-linear response of the calorimeter to particle 
energy, the loss of energy in uninstrumented regions 
of the detector, and the energy radiated outside of the 
jet cone. 



3. Event Selection 

Events are selected which have large $! T and exactly 
two jets whose invariant mass can be reconstructed. 
This signature is sensitive to both t vqq and vvqq de- 
cays because a charged lepton is not explicitly required 
in the final state. Due to the limited dijet mass resolu- 
tion, there is significant overlap between the W — > qq' 
and Z — > qq dijet mass peaks. Therefore we con- 
sider as signal the combination of three diboson sig- 
nals (WW, WZ, and ZZ), and search for diboson 
production in the dijet mass range 40 < Mjj < 160 
GeV/c 2 . 

Events are selected using a set of |5 T -based triggers 
with a variety of $! T and jet requirements. All of these 
triggers have benefited significantly from a calorimeter 
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trigger upgrade completed in 2007 [ID]. The majority 
(94%) of events satisfy the inclusive $ T trigger which 
requires $ T > 45 GeV. We require events to have 
iftrp > 60 GeV and exactly two jets with Et > 25 
GeV and I77I < 2.0, which ensures a trigger efficiency 
of 96% ± 2% on signal as measured in Z — > events. 
Additionally the fraction of the total jet energy which 
is in the electromagnetic calorimeter is required to be 
less than 90% to eliminate the possibility that elec- 
trons and photons are reconstructed as jets. Sev- 
eral criteria suppress contamination from non-collision 
backgrounds. Events are required to have at least one 
reconstructed vertex formed by charged particle tracks 
to remove cosmic-ray events. To reduce beam-related 
backgrounds, the electromagnetic fraction of the total 
event energy must be greater than 30%. Addition- 
ally, the arrival time of the jets as measured by the 
electromagnetic shower timing system must be 
consistent with the pp collision time. After all selec- 
tion criteria are made, non-collision backgrounds are 
a negligible contribution to the sample (fewer than 90 
events out of the total 44,910 selected events). 

In order to supp ress the multijet background de- 
scribed in Sec. |4.2[ we use a $ T resolution model to 
distinguish true ]fi T originating from neutrinos from 
fake Ifjrp arising from mismeasurement of jet energies. 
The Ifjrp significance is a dimensionless quantity based 
on the event topology, the energy resolution of the 
jets, and the soft unclustered particles. The $ T sig- 
nificance is calculated as follows: 



$ T significance = — log-^T 5 ), 
where V = / V l {x)dx if < 



x10 j 



or V = 



n 



Vi{x)dx if Vi > 0, 



with yi = $ T /(E i T cos (1) 



In Eq. [T] Vi(x) is the jet energy resolution function 
of the i-th jet, E % T is the energy of the i-th jet, and 
A(f>i is the azimuthal angle between the i-th jet and 
]f T . The jet energy resolution function V[x), defined 
as a ratio of the detector-level and hadron-lcvcl jet 
energy, is obtained from dijet pythia Monte Carlo 
[12] as a function of a jet's energy and pseudorapid- 
ity The function is validated in inclusive Z — > ee 
events. The ~$ T significance is typically low when $ T 
arises from mismeasurement. In addition to having a 
small significance, the ]fi T will often be aligned with a 
jet. Distributions of the $ T significance and azimuthal 
angle between the $ T and nearest jet (A^ff ) after 

all previously described selection criteria are shown 
in Figure [l] To reduce the multijet background, we 
select events with ]p T significance larger than 4 and 
Aftrf greater than 0.4 radians. 
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Figure 1: Data compared with the sum of the predicted 
electroweak (EWK) and multijet background (MJB) con- 
tributions for the $ T significance (top) and A0 J J1 (bot- 

torn) variables defined in Sec. (3] The expected signal is in- 
cluded here in the EWK contribution. The gray bands rep- 
resent the total systematic uncertainty described in Sec. [6] 
To reject events with fake $ T we require $ T significance 

> 4 and A0 j ^ > 0.4 radians. 
ifj T 



4. Data Modeling 



The diboson signals (WW, WZ, and ZZ) are sim- 
ulated by the pythia Monte Carlo generator. After 
selection, the most significant backgrounds to the di- 
boson signal are multijet production and electroweak 
(EWK) processes such as W{— * ^)+jets and Z(— * 
z^p)+jets. Less significant EWK backgrounds include 
Z(—> ££)+jets, top-quark pair production, and sin- 
gle top-quark production. The geometric and kine- 
matic acceptance for all electroweak processes are de- 
termined using a Monte Carlo calculation of the col- 
lision followed by a GEANT3-based simulation of the 
CDF II detector response [13] . Modeling of the back- 
grounds is described in more detail in the following 
sections. 
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4.1. Electroweak Background 

The electroweak (EWK) processes considered to be 
backgrounds to this measurement are the K+jet and 
top-quark production processes. The W(— > £D)+jets 
processes are simulated by the fixed-order matrix el- 
ement generator ALPGEN [T3] interfaced with pythia 
to simulate parton showering and fragmentation, the 
underlying event, and additional pp interactions in the 
event. The Z+jets and top-quark production are sim- 
ulated by pythia. The expected yields of all Monte 
Carlo-simulated background processes are normalized 
to SM cross sections calculated at next-to-leading or- 
der, although the overall normalization of this back- 
ground is allowed to float in the final fit to data. 



4.2. Multijet Background 

The multijet background (MJB) does not typically 
result in signatures of large $ T . However, when jet 
energy is not measured accurately an event may be 
reconstructed with large fake ]fi T and thus pass the 
event selection. Due to the large rate of MJB produc- 
tion, this is still a significant background. Much of 
the MJB is removed by the requirements on the Ap 
significance and AcfP^ variables described in Sec. 3 

The shape and normalization of the remaining MJr 
are determined from the data. 

To determine the MJB contribution, a region dom- 
inated by MJB is isolated in the data. A missing 
transverse momentum vector ft T , analogous to the 
calorimeter-based $! Tl is constructed from the vector 
sum of the transverse momenta of all particles mea- 
sured in the tracking system. For $J T arising from neu- 
trinos, the f!rp and $J T will usually be closely aligned 
and the azimuthal angle A(j)(]fi T ,$ T ) between them 
will be small. The f T is largely uncorrected to the $ T 
in events where jet energies have not been measured 
accurately; thus the MJB is expected to be the domi- 
nant contribution at large values of A(j)(]£ T ,fi T ). The 
distribution of A<f)(]ft T ,ft T ) for events satisfying the se- 
lection criteria is shown in Figure [2j The EWK Monte 
Carlo is normalized to the data in the peak region 
of A(f>(]p Tj f> T ) < 0.3, where the MJB contribution is 
negligible. Additional corrections to the Monte Carlo 
are made based on comparisons of Z — ► Monte 
Carlo to Z — > /i/i events in data. Because muons are 
minimum ionizing particles which deposit little energy 
in the calorimeter, from calorimeter activity Z — > /i/i 
events appear to have real Jf, T . This makes Z — -> /i/i 
an excellent "standard candle" for processes with ~$ T . 

The region of A(f>(]p T ,f T ) > 1.0 is dominated by 
MJB. Subtracting the corrected EWK Monte Carlo 
from the data in the A<j}{]p T ,tf T ) > 1.0 region gives an 
estimate of the remaining MJB in the signal sample. 
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Figure 2: Data compared with the sum of the predicted 
electroweak (EWK) contributions, including the expe cted 
signal, for the A<j>{$ T ,f T ) variable defined in Sec. 4.2 
The EWK contribution is determined from Monte Carlo, 
and is corrected by normalizing to the data in the peak 
region, A(f>($ T ,f T ) < 0.3. Events with A<t>{$ T ,f T ) > 1.0 
are used to model the multijet background. 



We now determine the shape of the MJB in any vari- 
able of interest using events in the A(f>(]£ T , f T ) > 1.0 
region, with the EWK contribution subtracted. The 
shapes of the derived MJB distributions are also ver- 
ified with a high-statistics dijet pythia Monte Carlo 
sample. The systematic uncertainties associated with 
this data-driven method of determining the MJB con- 
tribution come primarily from potential differences 
between the MJB contribution in the Acj)(]fi T ,f T ) < 
1.0 and Acf>(]p T ,f T ) > 1.0 regions. To estimate the 
size of these uncertainties, we compare distributions 
in these regions in the dijet Monte Carlo sample. The 
normalization must also be scaled up to account for 
MJB contamination in the region A<f)(]fi T ,'fi T ) < 1.0, 
and an uncertainty of 20% is applied to the total MJB 
yield. 



5. Analysis Technique 

The analysis proceeds via a fit of the signal and 
background contributions to the dijet mass Mjj dis- 
tribution. The signal extraction is performed with an 
unbinncd maximum likelihood fit using the RooFlT 
program |15j . Three Mjj template distributions are 
used in the fit. The first is the signal shape. This 
template is obtained from a Gaussian plus polyno- 
mial fit to the signal Monte Carlo Mjj distribution, 
which is a combination of the WW, WZ, and ZZ 
Mjj distributions weighted by the predicted SM cross 
sections. The mean and width of the Gaussian arc lin- 
early dependent on the jet energy scale (JES), which 
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is constrained to be in the range allowed by external 
measurements j9j. The second template is the EWK 
shape, a combination of the Mjj distributions for all 
backgrounds taken from Monte Carlo as described in 
Sec. 4.1 The third is the MJB template, which is 



determined by forming the Mjj distribution in events 
with A(j)(]fi T ,f T ) > 1.0 as described in Sec. 



4.2 



The 

MJB Mjj distribution is modeled by an exponential 
which is used as the template in the fit. Based on 
studies with the dijet PYTHIA Monte Carlo sample, a 
20% uncertainty is taken on the slope of the exponen- 
tial to account for uncertainties on the shape of the 
MJB contribution. 

In the fit to data, the JES, the MJB normalization, 
and the slope of the MJB exponential enter as Gaus- 
sian constraints, allowed to float within the predeter- 
mined uncertainties. The yield of signal and EWK 
background events are floating in the fit with no con- 
straints. 



6. Systematics 



Table I The systematic uncertainties described in Sec. [6] 
and their effect on the number of extracted signal events, 
the acceptance, and the cross section measurement. All 
systematics are added in quadrature. Shape uncertain- 
ties on the MJB background and JES are included in the 
statistical uncertainty returned by the fit. 





Systematic Source 


% Uncertainty 


Extraction 


EWK background shape 
Resolution 


7.7 
5.6 




Total extraction 


9.5 




JES 


8.0 




JER 


0.7 


Acceptance 


$ T resolution model 


1.0 




Trigger inefficiency 


2.2 




ISR/FSR 


2.5 




PDF 


2.0 




Total acceptance 


9.0 




Luminosity 


5.9 




Total cross section 


14.4 



We address separately two classes of systematic 
sources, those that affect the signal extraction and 
those that affect the signal acceptance in the cross 
section calculation. The signal extraction systematic 
uncertainties come from uncertainties on the shapes 
of the signal and background templates. These shape 
uncertainties include effects of the jet energy resolu- 
tion (JER), the jet energy scale (JES), and the un- 
certainties on the shapes of the MJB and EWK back- 
ground. The JES and the shape and normalization of 
the MJB are treated as nuisance parameters in the fit, 
with a Gaussian constraint to their expected values. 
Thus these uncertainties are already accounted for in 
the statistical uncertainty on the extraction. The sig- 
nal resolution uncertainty is determined by smearing 
each jet in the Monte Carlo by the expected uncer- 
tainty on the JER. 

The EWK background Mjj template is taken from 
Monte Carlo modeling of all non-MJB backgrounds. 
The primary EWK backgrounds include a gauge bo- 
son accompanied by jets. To determine a shape uncer- 
tainty on the EWK background, an alternative back- 
ground model is developed using 7+jets data. There 
are similarities between V+jets and 7+jets produc- 
tion, but the kinematics are not identical due pri- 
marily to the large mass difference between the W/Z 
bosons and the photon. Also, selection requirements 
applied to the 7+jets events are not identical to those 
used in the ^ T +jets sample. The W/Z decays involv- 
ing neutrinos leave a signature of $ T in the detector, 
while the photon's Et is measured in the calorimeter. 
For this reason, in 7+jets events we require the vec- 
tor sum of the photon Et and any $ T in the event to 
be greater than 60 GeV, treating this sum as anal- 
ogous to the frp in V+jets events. To correct for 



the differences in kinematics and selection, the 7+jets 
events are weighted by the ratio of the Mjj distri- 
butions in the EWK background Monte Carlo sam- 
ples to a 7+jets pythia sample. This method ac- 
counts for any production differences, while allowing 
common sources of uncertainty, such as detector ef- 
fects, parton distribution functions, and initial and 
final state radiation, to cancel out. Another consider- 
ation in the 7+jets data sample is the contribution of 
ry+V events, which will cause a peak in the Mjj distri- 
bution. We use a 7 + V pythia Monte Carlo sample 
to subtract this contribution. After these corrections, 
the difference between the 7+jets and default EWK 
background shapes is very small. We then use the 
adjusted 7+jets Mjj distribution to perform a signal 
extraction fit. The systematic uncertainty due to the 
shape of the EWK background is estimated as the 
difference in the final parameter values between the 
results obtained from this fit and the fit with the de- 
fault EWK background shape. This method accounts 
for a combined effect of JES, JER, and modeling of 
jets in the Monte Carlo. 

The dominant source of systematic uncertainty on 
the signal acceptance and the cross section measure- 
ment is the uncertainty associated with the jet energy 
scale. The JES affects several of the variables used 
in the event selection. The effect of the JES uncer- 
tainty is quantified by varying the jet energies in the 
signal Monte Carlo by the ±1<7 variations of the JES. 
Other sources of systematic uncertainty affecting the 
measured cross section are the JER, the $ T resolution 
model used in calculating the $ T significance, the trig- 
ger inefficiency calculated from Z — > events, the 
initial and final state radiation (ISR/FSR), and the 
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partem distribution functions (PDF). An additional 
uncertainty originating from the luminosity measure- 
ment is also taken into account [7] . A summary of the 
sources of systematic uncertainty is given in Table [T] 



7. Results 

The measured yields for signal and background from 
the Mjj fit to data are given in Table |TT| with the 
extracted number of signal events measured to be 
1516 ± 239(stat) ± 144(syst). The acceptances for the 
WW, WZ, and ZZ processes are estimated as 2.5%, 
2.6%, and 2.9% respectively from the Monte Carlo 
simulations. In the calculation of the combined cross 
section, we assume each signal contributes proportion- 
ally to its predicted SM cross section, which is 11.7 pb 
for WW, 3.6 pb for WZ, and 1.5 pb for ZZ. Since 
the sample utilizes a large number of trigger paths, the 
luminosity of the sample is best calculated by count- 
ing the number of Z — > events in the ^-triggered 
sample and comparing to the well understood muon- 
triggered data. This gives a total effective luminosity 
of 3.5 fb . With this information, the measured sig- 
nal events correspond to a cross section of o~{pp — > 
VV + X) = 18.0 ± 2.8(stat) ± 2.4(syst) ± l.l(lumi) 
pb, in agreement with a SM prediction of 16.8 ± 0.5 
pb obtained from the MCFM program |16j using the 
CTEQ6.1M PDFs [17]. 

Figure [3] shows the fit result and a comparison be- 
tween the extracted signal and the data after back- 
ground subtraction. With the data binned as in Fig- 
ure]^ we obtain a x 2 of 9.4 for 9 degrees of freedom, 
corresponding to a p- value of 40%. To assess the sig- 
nificance of the observed signal, we consider parameter 
variations for the sources of systematic uncertainty. 
Then the maximum likelihood value of a background- 
only fit is compared to the likelihood value returned 
by the full fit, and the difference is converted into 
significance numbers. This method finds that the ob- 
served signal corresponds to a significance of at least 
5.3 standard deviations from the background-only hy- 
pothesis. 

In summary, we use a $! T plus dijet final state 
to measure the combined WW + WZ + ZZ cross 
section in pp collisions at y/s = 1.96 TeV to be 
18.0 ± 2.8(stat) ± 2.4(syst) ± l.l(lumi) pb. This is 
consistent with the SM prediction, and is the first ob- 
servation in hadronic collisions of the clcctroweak pro- 
duction of vector boson pairs where one boson decays 
to a hadronic final state. 
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Figure 3: Top: Comparison between data and the results 
of a maximum likelihood fit, with the signal shown un- 
stacked. The gray band represents the systematic uncer- 
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boson signal (solid line) with the background-subtracted 
data (points). The dashed lines represent the ±1<t sta- 
tistical variations on the extracted signal, while the gray 
band again represents the systematic uncertainty due to 
the EWK background shape. 



Table II Value of parameters in the model used to fit the 
Mjj distribution. Uncertainties shown are statistical only. 
Of the 44,910 data events which pass the selection criteria, 
1516±239 signal events are extracted. The jet energy scale 
is also extracted from the fit and agrees with the default 
value of 1.0 as measured from calibrations [9]. 



Parameter 



Fitted Value 



Jet energy scale 0.985 ± 0.019 

Yield of EWK background events 36140 ± 1230 
Yield of MJB events 7249 ± 1130 

Yield of diboson candidates 1516 ± 239 
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